The numerical stability of leaping methods for stochastic simulation of chemically reacting systems.
نویسندگان
چکیده
Tau-leaping methods have recently been proposed for the acceleration of discrete stochastic simulation of chemically reacting systems. This paper considers the numerical stability of these methods. The concept of stochastic absolute stability is defined, discussed, and applied to the following leaping methods: the explicit tau, implicit tau, and trapezoidal tau.
منابع مشابه
Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method
We show how stiffness manifests itself in the simulation of chemical reactions at both the continuous-deterministic level and the discrete-stochastic level. Existing discrete stochastic simulation methods, such as the stochastic simulation algorithm and the ~explicit! tau-leaping method, are both exceedingly slow for such systems. We propose an implicit tau-leaping method that can take much lar...
متن کاملImproved Leap-Size Selection for Accelerated Stochastic Simulation
In numerically simulating the time evolution of a well-stirred chemically reacting system, the recently introduced “tau-leaping” procedure attempts to accelerate the exact stochastic simulation algorithm by using a special Poisson approximation to leap over sequences of non-critical reaction events. Presented here is an improved procedure for determining the maximum leap size for a specified de...
متن کاملImplicit Simulation Methods for Stochastic Chemical Kinetics
In biochemical systems some of the chemical species are present with only small numbers of molecules. In this situation discrete and stochastic simulation approaches are more relevant than continuous and deterministic ones. The fundamental Gillespie’s stochastic simulation algorithm (SSA) accounts for every reaction event, which occurs with a probability determined by the configuration of the s...
متن کاملConsistency and Stability of Tau-Leaping Schemes for Chemical Reaction Systems
We develop a theory of local errors for the explicit and implicit tau-leaping methods for simulating stochastic chemical systems, and we prove that these methods are first-order consistent. Our theory provides local error formulae that could serve as the basis for future stepsize control techniques. We prove that, for the special case of systems with linear propensity functions, both tau-leapin...
متن کاملStochastic simulation of chemical kinetics.
Stochastic chemical kinetics describes the time evolution of a well-stirred chemically reacting system in a way that takes into account the fact that molecules come in whole numbers and exhibit some degree of randomness in their dynamical behavior. Researchers are increasingly using this approach to chemical kinetics in the analysis of cellular systems in biology, where the small molecular popu...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- The Journal of chemical physics
دوره 121 24 شماره
صفحات -
تاریخ انتشار 2004